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Abstract:  An  investigation  was  made  of  the  heat-up,  ignition  and 

combustion  of  a  nitrate  ester  monopropellant  droplet  when 
suddenly  placed  in  a  high  temperature  gas.  At  moderate 
pressures  a  combustion  zone  became  noticeably  developed 
as  the  droplet  approached  its  wet  bulb  temperature  and  a 
steady  burning  period  followed.  At  high  pressures  reactive 
effects  became  noticeable  early  in  the  heat-up  period  and 
burning  was  not  steady.  In  this  regime  the  droplet  spent 
most  of  its  lifetime  in  the  preignition  period.  A  heat-up 
theory  of  ignition  is  developed  and  found  to  be  in  agree¬ 
ment  with  the  moderate  pressure  data.  Burning  rates  also 
reported  over  the  test  range. 
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? 
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k  -  Li  uid  Thermal  Conductivity  (BTU'FT  SC;CR- 

k  -  Gas  Thermal  Conductivity  (BTU/FT  SEC°R,' 
g 

L  -  Heat  of  Vaporization  (BTU/LBM/ 

1^  -  Minor  Axis  of  an  Ellipsoid  (FT) 

1^  -  Major  Axis  of  an  Ellipsoid  (FT) 

m  -  Mass  Burning  Rate  vLBM/SEC) 

Nu  -  Dimensionless  Nusselt  Number 
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X  -  Frequency  Factor  in  Arrhenius  Expression 
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a  *  Constant  *  Bi/2,  Dimensionless 

B  *  Constant  *  a  -  1  *  Bi/2  -  1,  Dimensionless 

X  »  L/C  T  *  Dimensionless 

p  oo 

a  *  Stefan-Bolt2man  Coefficient 

«  Dimensionless  Ratio  of  Flame  Radius  to  Drop  Radius 
0  =  Density  of  Droplet  (LSM^FT^) 
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-  Density  of  Gas  at  Flame  Surface  (LBM/FT  ) 
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CHAPTER  1 


INTRODUCTION 

1. 1  General  Statement  of  the  Problem 

In  recent  years  monopropellant  droplet  combustion  has 
received  an  increasing  amount  of  attention.  While  liquid  monopro- 
pellantp  have  not  been  found  to  be  promising  for  use  in  the  primary 
propulsion  units  of  aircraft  or  missilesj  they  may  be  adapted  to 
secondary  applications.  A  thorough  understanding  of  the  combustion 
processes  of  a  single  monopropellant  droplet  is  an  important  first 
step  in  comprehending  the  combust icn  process  of  these  fuels  in 
rocket  motors.  This  thesis  is  concerned  with  an  experimental  and 
analytical  investigation  of  the  heat-up  and  burning  rates  of  a 
particular  liquid  monoprcpel lant  when  suddenly  subjected  to  a  high 
temperature  environment. 

1.2  Previous  Related  Studies 

Regarding  droplet  combustion  studies  in  general.  Robayasi^ 
and  Nishivaki^  introduced  the  use  of  an  electrically  heated 
movable  furnace  in  their  studies  of  bipropellant  ignition.  The 
heated  furnace  was  moved  ever  a  droplet  supported  on  a  quartz 
fiber  to  provide  rapid  introduction  cf  the  droplet  into  the  high 
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temperature  environment.  Motion  pictures  were  taken  of  the  heat¬ 
up  and  ignition  of  the  droplet. 

12  17 

A  technique  similar  to  that  of  kobayasl  and  Nishivaki 

3 

was  applied  to  monopropellant  combustion  by  5arrere  for  a  variety 

3 

of  monopropellant6 .  Barrere  tested  in  a  nitrogen  environment  at 
atmosphere  pressure  with  ambient  temperatures  between  1000°?  and 
1600°F.  Droplets  with  diameters  between  .030  and  .080  inches  were 
tested.  However,  primary  emphasis  was  on  the  quasi- steady  state 
burning  period  of  the  droplet  and  no  ignition  measurements  were 
made. 

3 

Barrere' s  test  is  the  only  study  tc  date  in  which  the  mono¬ 
propellant  was  apparently  allowed  to  burn  in  an  inert  atmosphere 
as  is  the  case  in  an  engine  where  the  monopropellant  burns  in  its 
products  of  combustion.  Tarifa,  et  al,^  and  Rcsser^  were  able 
to  support  combustion  only  with  the  addition  of  oxygen  tc  the 
environment  thereby  creating  a  system  not  truly  represented  as 
monopropellant  combustion. 

An  early  attempt  at  theoretically  modeling  the  combustion 

of  a  motionless  monopropellant  droplet  was  undertaken  by  Lcrell 
14 

and  Wise  .  Their  model  considered  a  chemical  reaction  diffused 
throughout  the  gas  phase  with  the  fuel  vapor  decomposing  exother¬ 
mically  into  product  gases  at  a  rate  dependent  cn  the  temperature 
and  weight  fraction  of  the  fuel  vapcr.  Ey  taking  equal  values  of 
molecular  weights  and  by  assuming  a  Lewis  number  of  unity,  the 


problem  was  solved  by  Integrating  numerically  the  resulting 

differential  equations  of  the  process. 

26  28 

Later  works  of  Tarifa,  et  al  ,  and  Williams  considered 

diffusion  of  species  of  different  molecular  weights  and  Lewis 
28 

numbers.  Williams  presents  exact  results  for  his  model  by 

applying  the  assumption  of  a  thin  flame  zone  for  reactions  with 

sufficiently  large  activation  energies.  This  model,  first 

24 

attributed  to  Spalding  ,  allows  direct  simultaneous  solution  of 
the  system  equations.  However,  the  general  lack  of  experimental 
results  has  prevented  any  verification  or  revision  of  existing 
theories. 

1  - 3  Specific  Statement  of  the  Problem 

The  preceding  discussion  has  indicated  that  knowledge  of 
the  combustion  processes  of  liquid  mcnopropellants  i?,  at  best, 
very  limited. 

Concerning  the  existing  theories  of  moncpropellant  droplet 
burning,  it  would  be  of  interest  tc  determine  what  limitations, 
if  any,  they  present.  Without  this  knowledge,  the  accuracy  of 
burning  lifetimes  predicted  by  these  theories  is  subject  to 
question. 

Also,  it  is  of  importance  to  understand  the  heat-up  mechanism 
of  a  monopropollant .  Its  relationship  to  the  heat-up  time  and, 
ultimately,  to  the  total  droplet  lifetime  must  be  a  primary  consid¬ 


eration  in  rocket  engine  design. 
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The  technique  employed  by  Barrere  offers  a  means  of 
studying  the  droplet  combustion  process.  However,  by  enclosing 
the  apparatus  in  a  pressure  vessel  equipped  with  a  means  of  evacua¬ 
tion,  a  pure  inert  atmosphere  may  be  maintained  at  various  total 
pressures . 

To  summarize,  the  specific  objectives  of  this  thesis  are 
as  follows: 

1.  To  formulate  a  theory  for  the  heat-up  and  ignition  of 
a  motionless  monopropellant  droplet,  and  check  ..  experimentally 
over  a  range  of  ambient  temperatures  and  pressures. 

2.  To  apply  a  suitable  existing  mor.opropel lant  burning 
theory  to  the  experimental  results,  and  thereby  determine  its 
limitations. 
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CHAPTER  2 


EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

2 . 1  Test  Facility 

The  primary  requirement  of  the  experimental  apparatus  was 
to  provide  a  means  of  rapidly  placing  a  droplet  into  a  high 
temperature  inert  atmosphere  at  various  total  gas  pressures.  A 
means  of  photographically  recording  the  droplet's  behavior  and 
a  means  of  measuring  the  droplet's  temperature  were  also  required. 

A  photograph  of  the  experimental  apparatus  used  for  this 
study  of  the  ignition  and  combustion  of  monopropellant  droplets 
is  shown  in  Figure  2.1.  Figure  2.2  is  a  schematic  diagram  of  the 
internal  chamber  assembly. 

The  upper  part  of  the  apparacus  was  enclosed  with  a  pressure 
vessel  which  allowed  operation  up  to  3,000  PSI.  A  hand-operated 
winch  was  used  to  raise  and  lower  the  upper  cover  and,  when  the 
cover  was  bolted  to  the  base,  the  chamber  was  gas  tight.  By  using 
a  vacuum  pump  and  a  removable  hose  attached  to  a  pressure  sealed 
connector,  the  chamber  was  evacuated  and  then  pressurized  with  the 
inert  gas.  For  operating  pressures  below  one  atmosphere,  the 
pressure  within  the  chamber  was  measured  with  a  differential 
mercu- y  manometer.  Measurement  of  chamber  pressures  between  i  and 
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5  atmospheres  was  made  with  an  Ashcroft  laboratory  gauge  with 
1/2  PSI  subdivisions.  For  tests  above  5  atmospheres,  an  Ashcroft 
gauge  with  50  PSI  subdivisions  was  used.  Conmercially  pure 
nitrogen  was  used  to  pressurize  the  chamber.  The  monopropellant 
used  in  the  testing  program  was  propylene  glycol  dinitrate. 

The  droplet  was  supported  on  either  a  quartz  filament  or 
the  bead  of  a  thermocouple  junction.  A  sketch  of  the  mount  assembly 
is  shown  in  Figure  2.3.  The  thermocouple  was  constructed  from 
chrome 1-alumel  wire  having  an  outside  diameter  of  .003  inches. 

The  somewhat  enlarged  thermocouple  junction  had  an  average  diameter 
of  between  .008  and  .016  inches. 

When  only  one  thermocouple  was  used,  it‘s  output  was  fed 
directly  into  the  vertical  terminals  of  a  Model  130B  Hewlett- 
Packard  oscilloscope.  A  Dumont,  Type  2614,  oscilloscope  camera 
employing  Type  47  Polaroid  film  was  used  to  record  the  trace.  In 
the  dual  thermocouple  arrangement,  the  thermocouple  outputs  were 
fed  into  the  terminals  of  two  galvanometers  in  a  Consolidated 
Electrodynamics  Corporation  Model  5-124  recording  oscillograph. 

The  quartz  filament-single  thermocouple  arrangement  was  used 
primarily  in  the  larger  droplet  sizes  (.060  -  .080  in.;  because 
the  small  thermocouple  junction  could  not  support  the  large  droplets. 

The  furnace  consisted  of  a  heating  coil  wound  between  an 
outer  stainless  steel  cover  and  an  inner  ceramic  core.  The  furnace 
was  cylindrical  in  shape  with  an  inside  diameter  of  2  in.  and  an 
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outside  diameter  of  3-1/2  in.  The  over-all  height  of  the  furnace 
was  c-l/2  in.  with  an  actively  heated  inside  length  of  5  in. 

After  the  furnace  was  heated  to  the  desired  temperature, 
a  solenoid  was  energized  releasing  it  from  its  initial  position. 

The  furnace  slid  down  guide  rods  to  the  droplet  location  where  it 
was  locked  in  place  by  two  latches,  thus  providing  rapid  immersion 
of  the  droplet  into  the  high  temperature  media.  The  time  required 
for  the  furnace  bottom  to  travel  from  the  droplet  location  to  the 
latches  was  about  50  ms.  The  shock  of  the  fall  was  absorbed  by 
two  pieces  of  Resilite  mounted  on  a  false  base.  Resiilte  is  an 
energy  absorber  sometimes  used  in  athletic  equipment  to  reduce 
shock.  The  Resilite  shock  absorbers  combined  with  the  latching 
mechanism  eliminated  extraneous  motion  of  the  gas  around  the  drop¬ 
let  during  the  test  period. 

The  bottom  of  the  furnace  was  open,  while  the  top  was 
closed  with  a  quartz  window.  This  window  allowed  the  passage  of 
background  illumination  for  photographic  purposes,  while  preventing 
the  escape  of  high  temperature  gas  f  cm  the  furnace  irtc  the 
chamber. 

The  gas  temperature  within  the  furnace  was  measured  with  a 
1/8  in.  O.D.  sheathed  chrcrael  alumel  thermocouple  Tne  thermo¬ 
couple  junction  was  located  in  a  position  that  would  correspond 
to  the  droplet  position  when  the  furnace  was  locked  in  place. 

This  was  about  2-1/2  in.  from  the  base  of  the  furnace.  The 


11 


temperature  indicated  by  this  thermocouple  just  before  each  test 
was  taxen  as  the  test  temperature.  The  thermocouple  was  allowed 
to  stabilize  at  a  constant  temperature  before  each  test  to  insure 
a  uniform  furnace  temperature. 

The  droplet  was  photographed  through  a  quartz  window  in 
the  base  of  the  chamber  using  a  16mm  Fastair  missile  camera. 

Kodak  Tri-X  negative  film  was  used  in  the  camera  at  a  speed  of 
approximately  SO  frames  per  second.  The  background  light  was 
a  small  incandescent  lamp.  This  illumination  passed  through  the 
quartz  window  in  the  top  of  the  furnace  giving  a  shadowgraph  of 
the  droplet.  A  neon  timing  light  in  the  camera  allowed  synchroni¬ 
zation  of  the  motion  picture  film  and  the  oscilloscope  trace.  As 
the  furnace  reached  the  droplet  location,  it  closed  two  switches 
which  simultaneously  completed  the  circuits  of  the  timing  marker 
and  the  oscilloscope  trace,  thereby  indicating  the  beginning  of 
the  heat-up  process. 

2.2  Operation  of  the  Apparatus 

The  ovet-all  operation  of  the  apparatus  is  illustrated  in 
Figure  2.4.  Conditions  within  the  chamber  at  the  beginning  of 
a  test  are  represented  in  the  first  diagram  of  Figure  2.4.  At 
this  point  the  furnace  temperature  has  stabilized,  the  droplet 
has  been  placed  on  its  support,  and  the  chamber  hes  been  closed 
and  pressurized.  The  background  light  has  been  turned  on  and  the 
camera  is  off. 
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To  begin  operation  the  oscilloscope  sweep  was  triggered, 
the  camera  was  started,  and  the  furnace  release  solenoid  was  actu¬ 
ated  allowing  che  furnace  to  fail.  In  the  second  diagram  of 
Figure  2.4,  the  bottom  edge  of  the  furnace  is  just  passing  the 
droplet  location.  At  this  point  the  furnace  closed  two  switches 
which  simultaneously  completed  the  oscilloscope  trace  circuit  and 
the  timing  light  circuit  on  the  camera 

The  third  diagram  of  Figure  2.4  shews  the  furnace  in  its 
final  locked  position  at  the  droplet  location. 


CHAPTER  3 


IGNITION 


3 . 1  Preliminary  Results 

The  exact  nature  of  the  ignition  process  of  a  particular 
fuel  varies  with  the  chemical  composition,  environment,  etc  of 
the  fuel  being  investigated.  Consequently,  preliminary  experi¬ 
mental  tests  seemed  advisable  to  assist  in  the  selection  of  the 
ignition  model.  Typical  results  of  these  and  later  tests  are 
shown  in  Figures  3.1  and  3.2. 

Figures  3.1a  and  3.1b  which  present  da~a  at  1  atm  total 
pressure  over  the  experimental  temperature  range  demonstrate  a 
correlation  between  the  beginning  of  diameter  change  and  an  inflec¬ 
tion  point  in  the  liquid  phase  thermocouple  record.  The  diameter 
change  is  the  point  in  the  diameter  versus  time  plct  where  the 
diameter  begins  to  decrease  rapidly.  It  occurs  very  slightly 
before,  or  about  the  same  time  as  the  inflection  point.  A  gas 
phase  thermocouple  was  then  introduced  to  indicate  when  a  flame 
appeared  in  the  gas  phase,  since  the  flame  was  not  visible  at 
low  pressures.  It  is  seen  that  the  inflection  in  the  gas  phase 
thermocouple,  which  indicates  the  presence  cf  a  flame,  exhibits 
the  same  correlation  with  the  diameter  change  that  the 
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thermocouple  had,  as  evidenced  in  Figure  3,2b.  That  is,  the 
liquid  and  gas  traces  inflect  at  the  same  time,  or  very  slightly 
after,  the  diameter  change  occurred.  This  correlation  extends 
over  the  entire  test  range.  Figure  3.2a,  which  presents  data 
at  40  atmospheres  pressure,  shows  excellent  agreement  between 
the  inflections  of  the  two  curves  and  the  appearance  of  the 
luminous  flame.  The  diameter  of  the  droplet  during  the  pre¬ 
ignition  period  could  not  be  measured  accurately  at  high  pressures, 
due  to  the  scattering  of  the  background  light  rays  by  natural 
convection  patterns  on  the  furnace  wtudov. 

On  the  basis  of  the  preceding  Intones,  ion,  a  criterion 
for  the  ignition  time  was  established.  Each  acceptable  data 
point  had  to  have  at  least  two  separate  indications  oj  ignition 
that  occurred  at  very  nearly  the  same  time.  Gas  or  liquid  thermo¬ 
couple  inflection  and  diamet^i.  change  or,  at  high  pressures,  gas 
or  liquid  thermocouple  inflection  and  luminous  flame  were  the 
combinations  used. 

Figures  3.3  and  3.4  show  sample  photographs  of  the  ignition 
and  combustion  processes.  They  illustrate  the  extremes  of  the 
problem  of  ignition  identification.  In  Figure  3.3,  where  the 
visible  flame  was  used  as  the  ignition  indicator,  the  film  record 
between  .31  and  .85  seconds  was  eliminated  since  no  significant 
changes  in  the  droplet  were  observed  during  that  period.  Note 
how  the  natural  convection  patterns  on  the  furnace  window  in 


SAMPLE  PHOTOGRAPH  -  VISIBLE  FlAME 


PHOTOGRAPH  -  FLAME  NOT  VISIBLE 
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Figure  3.3,  at  40  atmospheres  total  pressure,  tend  to  obscure 
the  droplet  during  the  pre- ignition  period.  No  such  problem 
was  encountered  in  the  1  atmosphere  test  of  Figure  3.4.  The 
droplet  regained  essentially  unchanged  between  .17  and  .88  seconds, 
so  this  film  was  excluded.  Shortly  after  .88  seconds  the  droplet 
begins  to  steadily  decrease  in  size  until  it  is  almost  completely 

consumed  at  1.03  seconds.  The  empty  probe  is  all  that  remains 

3  7  12 

at  1.40  seconds.  As  noted  by  other  experimenters,  *  *  the 
droplet  is  not  exactly  spherical  during  the  test  period.  To 
retain  a  degree,  of  uniformity  with  previous  workers  in  the  field, 
the  droplet  was  assumed  to  nave  the  shape  of  an  ellipsoid.  The 
diasieter  of  the  droplet  was  then  defined  to  be  the  diameter  of  a 
sphere  having  a  volume  equal  to  the  volume  of  the  ellipsoid. 

The  equivalent  sphere’s  disaster  is  related  to  the  minor 
and  major  axes,  and  l 2,  of  the  ellipsoid  by  the  equation 


D 


(ii2  V 


1/3 


Upon  examining  the  behavior  of  the  droplet  in  Figure  3.1,  a 
gradual  increase  in  the  rate  of  decrease  of  the  droplet  diameter 
siey  be  observed  after  ignition  has  occurred.  Two  possible  causes 
of  this  are  the  increased  strength  of  the  flame  as  the  burning 
process  progresses  and  the  increased  mass  transfer  rates  at  high 
vapor  pressures,  irrespective  of  flame  location. 
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The  liquid  temperatures  of  che  droplet  of  Figure  3-1 
apparently  begin  to  heat  up  to  a  wet  bulb  temperature  consistent 
with  the  temperature  of  the  environment.  However,  when  ignition 
occurs  the  droplet  sees  the  higher  temperature  of  the  flame  zone 
and  its  temperature  then  rises  to  an  equilibrium  temperature  in 
the  presence  of  the  flame. 

3.2  Theoretical  Considerations 

On  the  basis  of  these  preliminary  tests,  two  theoretical 
models  were  developed  to  describe  the  pre- ignition  behavior  of 
a  motionless  droplet  instantly  surrounded  by  a  high  temperature 
inert  atmosphere.  The  following  assumptions  were  made  in  the 
analysis : 

1.  The  droplet  is  composed  of  a  single  chemical  species. 

2.  In  the  first  model  the  droplet  temperature  is  assumed 
uniform,  but  varying,  while  the  second  model  considers 
transient  temperature  gradients  in  the  droplet. 

3.  Properties  are  constant  in  both  the  liquid  and  gas  phases. 

4.  No  surface  evaporation  occurs,  and  the  droplet  radius 
remains  constant. 

5.  The  heat  transfer  and  diffusion  processes  are  spheri¬ 
cally  symmetrical. 

6.  The  influence  of  mass  transfer  on  the  heat  transfer 
characteristics  of  the  boundary  layer  around  the 
droplet  is  neglected. 
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Since  ignition  occurred  when  the  droplet  temperature 
approached  the  wet  bulb  tempera ture,  which  is  relatively 
near  to  the  boiling  point  at  the  total  pressure  cf  the 
test,  it  appeared  reasonable  to  take  the  time  required 
for  a  droplet  to  heat  up  to  the  boiling  point  as  an 
estimation  of  ignition  time. 


Considering  first  the  uniform  temperature,  or  infinite 
conductivity,  model,  the  energy  equation  is 


oB(T4  -  T4)A  i-  hA;T  -  T)  =  V7 

«>  '  OC  pdc 


4  4 

assuming  T  »  T  ,  Equation  3.1  may  be  integrated  directly  to 
yield 
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The  Nusselt  number  is  determined  by  employing  "he  results  of 


Ranz  and  Marshall  for  natural  convection  heat  transfer  to  drcps. 
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Nu 


2.0  -i-  0.6  G, 


0.25  „  0.33 


Consider  now  the  finite  conductivity  case  with  temperature 
gradients  existing  in  the  droplet  In  this  analysis  the  complete 
energy  equation  must  be  considered  within  the  droplet.  It  is 
assumed  in  the  formulation  that  radiation  is  absorbed  uniformly 
throughout  the  droplet.  This  assumption  seems  realistic  for 
droplets  of  the  sizes  studied  in  view  of  the  theoretical  considera- 

9 

tions  by  Hottel  et  al  on  droplet  radiation  absorption  characteris- 
4  4 

tics.  Since  T  »  T  ,  the  surface  radiation  flux  is  taken  as 

CO 

4 

aBT  and  the  radiation  absorbed  per  unit  volume  is 
00 

3cBT4 

OO 

Q  a  - —  3  4 

r 

o 

The  energy  equation  within  the  droplet  for  the  finite  conductivity 
case  becomes: 


oC 
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P  dt 
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dt 


3.5 


The  boundary  conditions  available  for  the  solution  of  the 
energy  equation  are  the  surface  convection  condition,  and  the 
existence  of  a  finite  solution  at  the  droplet  center. 


24 


dT 

dr 


h(7  -  1), 

nr  -  * 


r  =  r„ 
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T  is  finite,  r  -  0 

The  equations  were  first  put  in  dimensionless  form  by  the 
following  substitutions; 

I  -  T 


T  -  7 
*  c 


1  =  r/r 


kt 


k-C  r 
P  o 


With  these  substitutions.  Equation  3.^  becomes 


£6 
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3.7 


Where 


Qr, 


k(T  -  T 

<r  o 
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and  the  ncn-dimensionali2ed  boundary  conditions  are: 


S’! 


n 


1 


d  is  finite,  q  =  0 


The  Laplace  transform  method  was  used  to  solve  Equation  3.7 
with  the  boundary  conditions  of  3.8.  The  complete  details  of  the 
solution  may  be  found  in  Appendix  A.  Tne  solution  fi  he  tempera¬ 
ture  distribution  in  the  sphere  as  a  function  of  time  with  connec¬ 
tion  and  radiation  heating  is 


2 


This  solution  was  evaluated  on  the  IBM  7074  computer  using 

23 

the  first  five  roots  of  Tan  Y  .  as  tabulated  by  Schneider  ,  for 

n 

various  values  of  Bi  and  \|r. 


2' 


The  radiation  absorption  coefficient  B,  for  both  models 

was  taken  to  be  0.5.  On  the  basis  of  the  experimental  results 

7  9 

of  Faeth  ,  and  the  theoretical  work  of  Hottel,  et  al  this  value 

seems  realistic  for  the  droplet  sizes  studied.  Property  values 

employed  in  the  calculations  may  be  found  in  Appendix  B. 

The  effect  of  including  radiation  in  the  finite  conductivity 

analysis  is  shown  in  Figure  3.5.  At  the  surface  of  the  droplet, 

the  temperatures  predicted  by  the  finite  conductivity  model  with 

no  radiation  are  15  to  20  per  cent  below  the  temperatures  of 

the  model  that  included  radiation.  Since  the  surface  temperature 

is  of  primary  importance  in  determining  the  theoretical  ignition 

delay  time,  it  is  apparent  that  radiation  should  be  included. 

Figure  3.5  also  indicates  that  significant  temperature  gradients 

exist  in  the  droplet.  The  importance  of  this  is  recognized  in 

relation  to  the  liquid  phase  thermocouple  readings.  Since  the 

thermocouple  is  located  at  some  point  in  the  droplet  between 

the  center  and  the  outer  edge,  it  will  indicate  a  temperature 

somewhat  below  the  actual  surface  temperature. 

Finally,  the  temperature  at  which  ignition  occurred  was 

assumed  to  be  the  boiling  tempera t  <re  of  the  fuel  droplet.  The 

boiling  temperatures  were  determined  by  employing  the  vapor 

4  6 

pressure  data  of  Beilstein  and  de  C  Crater  for  propylene  glycol 
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dinitrate  between  10  and  10  atms.  Reid  and  Sherwood  indicate 

that  a  plot  of  In  P  versus  1/T  usually  gives  a  reasonably  straight 


DIMENSIONLESS  TEMPERATURE  0  =  (T-To>/{ 


line.  By  fitting  a  straight  line  through  rhe  data  available,  and 
extending  it.  at  constant  slope,  to  higher  pressures,  the  bciiing 
temperatures  vere  obtained  at  those  pressures  The  straight  line 
assumption  of  Figure  3  6  seems  reasonable  m  view  of  Figure  3.7. 

All  three  compounds  plotted  using  Lange's  data  exhibited  the 
straight  line  characteristic  up  cc  the  critical  point. 

3 . 3  Experimental  Results  and  Discussion 

In  order  to  evaluate  th-  suitability  cf  the  prcpcsed 
ignition  models,  further  tests  vere  conducted.  The  results  cf  the 
experimental  program  are  presented  in  Figures  3  8  through  3.12 

Figure  3  8  presented  the  variation  of  ignition  delay  time 
with  initial  droplet  diameter  at  a  total  pressure  of  1  atmcspr.ere 
and  an  ambient  temperature  of  !^60CR  Although  toe  infinite  con¬ 
ductivity  model  is  in  rcugb  agreement  with  the  data,  the  finre 
conductivity  model  gives  a  superior  fit  of  the  r°nds  of  t^is 
figure  However,  due  tc  the  approximate  nature  of  the  node-; 
and  rhe  inaccuracies  due  tc  series  t-uncaticn  cf  the  fuutr 
conductivity  model  for  small  -im^s.  it  sterns  suffiCLenr  to  be;- 
the  comparison  c~,  toe  infinite  conductivity  model 

The  variation  or  ignition  deiev  ■  ime  vi'*<  ambient  temperature 
for  a  given  droplet  siee  is  presented  in  Fig-re  3  9  Again  the 
analysis  foli-'ws  the  general  'rend  of  the  experimental  results 
reasonably  v-rii.  The  change  m  igni'i  n  d-lav  is  ^vpisnsbli  bv 
faster  heat-up  rates  aT  higner  ambient  Temperatures  This  mdica't.s 
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that  ambient  temperature  exerts  primary  influence  on  heat-up  rates 

at  this  pressure,  rather  than  pre-ignition  chemical  kinetics  as  is 

the  case  for  bipropellant  ignition,^ 

Finally,  the  effect  of  various  total  pressures  on  the 

ignition  delay  time  is  presented  in  Figure  3.10  and  3.11.  As 

the  ambient  temperature  increases,  the  data  follows  the  analysis 

more  closely  at  low  pressures.  However,  above  1  atm.  total 

pressure  the  analysis  appears  to  break  down  soroewhar.  Considering 

Figure  3.11,  it  may  be  noted  that  in  all  three  cases  the  temperature 

indicated  by  the  liquid  thermocouple  falls  between  the  center  and 

surface  temperatures  of  the  finite  conductivity  model.  From  the 

0.1  atm  and  1  atm  data,  it  appears  reasonable  to  assume  that 

ignition  occurred  at  a  temperature  reasonably  close  to  the  wet 

bulb,  or  boiling  temperatures  predicted  by  Figure  3.7  However- 

the  ignition  temperature  of  the  11  atm  test  appears  to  be  far  from 

the  1060°R  temperature  estimated  from  Figure  3.7.  Figure  3.12 

shows  that,  in  the  lower  temperature  range,  the  maximum  measured 

liquid  temperature  compares  reasonably  ~ith  the  boiling  temperature 

determined  from  the  vapor  pressure  curve.  However,  above  1 

atmosphere  there  is  a  marked  departure  from  the  curve.  It  is 

possible  that  ignition  did  not  occur  at  the  boiling  temperature, 

but.  more  accurately,  at  the  wet  bulb  temperature  for  a  given 

pressure.  While  the  wet  bulb  and  boiling  temperatures  roughly 

correspond  at  low  pressures,  this  is  not  necessarily  the  case  at 
high  pressures. 


IN  DELAY  AT  VARIOUS  PRESSURES 
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Also,  it  is  important  to  note  that  in  none  of  the  data  at 
10  atmospheres  and  above  was  a  steady  state  burning  condition 
attained.  Apparently  the  flame  becomes  strong  enough  to  increase 
reaction  rates  even  when  fuel  concentrations  are  below  those 
which  would  be  attained  at  the  wet  bulb  condition.  Thus,  the 
burning  lifetimes  of  the  droplets  in  this  pressure  range  were 
very  small  in  comparison  to  the  heat-up  time,  and  the  combustion 
process  did  not  occur  at  a  uniform  temperature. 


CHAPTER  4 


COMBUSTION 

4.1  Theoretical  Considerations 

The  burning  rate,  or  mass  of  material  per  second  leaving 
a  droplet,  can  be  written  from  a  simple  mass  balance,  as 

d  A  *3  \ 

■  "  31  (j 5  <V  ; 


assuming  constant  density 


-  ■  f 


d-D-li 

dt 


4.2 


For  bipropellant  diffusion  flames  dD“/dt  =  -C  (a  constant'-. 

2  3 

This  D  correlation  was  used  by  Barrere  and  Moutet  in  cneir 

study  of  other  monopropellants,  so  in  order  tc  compare  results 

this  data  was  correlated  in  a  similar  manner..  It  is  recognized 

26  2 

from  the  results  of  Williams  that  the  D  relation  is  not  strictly 

valid  for  monopropellant  flames  However,  due  to  the  limited 

experimental  diameter  range  it  was  difficult  to  distinguish  whether 

dD/dt  *  const.,  or  dD*/dt  =  const  ,  or  whether  some  intermediate 

2 

case  existed,.  Using  the  D  correlation,  Equa! ion  4.2  becomes 
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m  *  ■  |  o  D  C  4.3 

where  C  is  known  as  the  burning  rate  constant  during  the  quasi¬ 
steady  state  burning  period.  The  value  of  C  for  a  particular 

2 

combustion  process  is  obtained  from  the  slope  of  the  D  versus 
time  plot  during  the  burning  period.  However,  since  the  droplet 
diameter,  D„  is  constantly  decreasing  during  the  burning  period. 
Equation  4.3  is  valid  at  any  instant  in  the  burning  period. 
Therefore,  to  obtain  a  moss  burning  rate  from  the  data  available, 
it  wa6  necessary  to  define  E  as  the  average  droplet  diameter 
during  the  quasi-steady  state  burning  period 

m  *  -  |  p  5  C  4.4 

The  problem  of  predicting  the  mass  burning  rates  of  menf:- 
propeliant  droplets  has  been  iir-fc sfigated  analytically  by  several 
authors  .  since  the  aialysis  presented  by  Williams  , 

for  motionless  droplets,  involves  measurements  and  properties 
readily  available  and  can  be  applied  to  nen-adi ibatic  burning 

it  appeared  most  useful  for  this  study.  The  simplified  model 

26  22 
presented  by  Williams  ,  and  first  attributed  to  Spalding*" 

i /eludes  the  following  assumptions. 
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hi: 


1.  The  activation  energy  is  sufficiently  large  sc  that 
all  heat  release  will  occur  in  a  thin  shell  at  some 
distance  from  the  droplet  surface. 

2.  The  droplet  temperature  is  constant 

3.  The  flame  is  spherically  symmetrical  and  natural 
convection  is  neglected 

4.  Thermal  diffusion  is  neglected 

5.  Quasi-steady  burning  is  considered. 

6.  All  gases  obey  the  ideal  gas  equation  of  state,  and 
the  total  pressure  is  constant 

The  equations  presented  bv  Williams^”  are 
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The  unknowns  in  the  above  equation  are  7  ,  m. 


and  p  v  . 
c  c 


However*  from  a  strictly  phenomenological  approach  tc  the  co^-sticti 
process,  assuming  a  first  order  reaction  with  one  step  kinetics 
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This  is  the  additional  relationship  available  to  solve  the 
system  of  equations. 

18  19 

The  experimental  and  theoretical  works  of  Phillips  ’  , 

Levy^  and  others**^  indicate  that  the  gaseous  decomposition  of 

f  E  \ 

organic  nitrates  obeys  the  Arrhenius  expression  K  *  X  EXP  ( 


^"r  T  / 
\  o  c 


The  basic  mechanism  of  decomposition  is  the  breaking  cf  the 

O^NO^  bond.  For  propylene  glycol  dinitrate,  as  reported  by 
18 

Phillips  ,  this  process  requires  an  activation  energy  of 
iced  1 

37.4  Thus,  although  the  fuel  obeys  the  first  order  law 

consistent  with  the  analysis,  its  activation  energy  may  not  be 
sufficiently  high  to  validate  the  thin  flame  assumption. 


4.2  Experimental  Results  and  Discussion 

Figures  4,1  and  4.2  present  the  mass  burning  rates  obtained 
by  correlating  the  test  data  with  Equation  4.3.  The  curves  fitted 
through  the  points  of  Figures  41  and  4  2  indicate  the  general 
trends  of  the  results  and  do  not  represent  the  predictions  of  an 
analytical  model.  1'ne  data  of  Figure  4  1  was  found  to  be  of  the 
same  numerical  level  as  the  data  cf  Barrere  and  Moutet^  for  ethyl 
nitrate  and  propyl  nitrate  in  the  same  temperature  range.  As  may 
be  noted  in  Figure  4  2,  limited  data  was  obtained  at  higher 
pressures.  This  was  due  to  the  fact  that  the  luminous  flame  present 
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in  that  region  prevented  measurement  of  the  droplet  during  the 

burning  period.  Consequently,  Equation  4.3  could  not  be  evaluated. 

It  should  also  be  noted  that  quasi-steady  state  burning  was  not 

attained  in  the  11  atm  results  cf  Figure  4  2,  as  the  droplet 

temperature  was  constantly  increasing  during  the  burning  process. 

Had  a  new  wet  bulb  temperature  been  maintained  during  the  combustion 

process,  the  mass  transfer  rates  would  probably  be  higher  The 

general  trend  of  the  results  of  Figures  4.1  and  4  2  is  for  the 

burning  rate  to  increase  with  increasing  ambient  temperature  and 

pressure.  This  is  consistent  with  Equation  4.8  since  Pcvc  may  be 

thought  of  as  a  burning  rate  per  unit  area  However,  it  does  not 

seem  reasonable  that  the  rapid  change  of  the  burning  rates  cf 

Figure  4.1  can  be  explained  by  ambient  temperature  changes  alone 

28 

To  compare  the  analysis  presented  by  Williams  ,  Equations 
4.5,  4.6  and  4.7  were  solved  simultaneously  on  the  IBM  7074  com¬ 
puter,  The  experimental  mass  burning  rate?  were  used  for  m  and 

the  three  unknowns  then  were  T  .  j  ,  and  a  v  If  the  combustion 

c  c  c  c 

process  followed  the  analysis  closely,  it  shculd  be  possible  to 

determine  the  activation  energy  from  a  log^v^  versus  1/T^  plot 

and  compare  it  with  the  previously  measured  activation  energy 
kca 1 

of  37.4  --T-  .  This  was  done  m  Figure  4.3  for  the  variable 

mole 

temperature  data  of  Figure  4.2,  using  varicus  vOje_  for  the 

thermal  conductivity.  It  was  found  that  i*.  =  nee-- sary  tc  use 

an  unrealistically  high  value  cf  the  nduct to  get 

^  i 

an  activation  energy  close  to  37. T  •  :  -1  possible  cau?=s 

BCi» 
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FIGURE  4.3 


ACTIVATION  ENERGY  ASSUMPTION 
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of  this  deviation  may  be  pointed  out.  First,  as  previously 
noted,  the  thin  flame  model  of  Williams  may  be  in  error.  The 
more  general  case  of  a  flame  of  finite  tnickness  would  result 
in  a  lower  flame  temperature  and  lower  mass  transfer  rates. 

Also,  it  is  important  to  note  that  adiabatic  burning  was  not 

g 

attained.  Faeth's  analysis,  which  considers  non- ad iabatic 
burning,  shows  that  the  non-adiabat ic  boundary  influences  the 
burning  rate  of  the  droplet  in  several  ways.  First,  the  heat 
transported  to  the  surroundings  is  not  available  to  vaporize  the 
fuel,  and  secondly,  the  reduced  fiame  temperature  decreases 
reaction  rates  which  causes  the  reaction  zone  to  move  away  from 
the  droplet.  Consequently,  the  temperature  and  concentration 
gradients  at  the  droplet  surface  are  decreased  and  a  corresponding 

O 

decrease  in  transfer  rates  results.  Finally,  Faeth's  analysis 
showed  that  the  flame  was  particularly  vulnerable  to  convection 
effects  for  non-adiabatic  burning.  Natural  convection  could 
cause  larger  heat  losses  from  the  flame  ani  tnus  explain  the 
abnormally  high  effective  thermal  conductivities  necessary  to 
get  reasonable  activation  energies.  These  results  aiso  help  to 
explain  the  data  of  Figure  4.1.  As  the  ambient  temperature  in¬ 
creased,  more  heat  was  supplied  to  the  droplet  Convection 
effects  on  the  reaction  zone  were  reduced,  and  the  reaction  zone 

moved  closer  to  the  droplet.  This  behavior  was  also  predicted 
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from  the  results  of  Williams  model,  since  j  decreased  with 

c 

increasing  ambient  temperature 


CHAPTER  5 


SUMMARY 

An  experimental  study  of  the  ignition  and  combustion  of  a 
liquid  monopropellant  in  an  inert  atmosphere  was  undertaken.  The 
primary  objectives  of  the  study  were: 

1.  To  determine  the  effect  of  various  ambient  temperatures 
and  pressures  on  the  ignition  delay  time  and  burning  rate  of  the 
monopr ope 1 la  nt . 

2.  To  propose  a  realistic  ignition  model  and  compare  it 
with  experimental  results. 

3.  To  check  the  validity  of  a  suitable  monopropellant 
burning  theory  when  applied  to  the  experimental  results. 

An  experimental  apparatus  was  employed  which  allowed  testing 
over  a  range  of  ambient  temperatures  and  pressures. 

On  the  basis  of  preliminary  tests,  two  heat-up  models  were 
proposed.  The  first  model  assumed  the  temperature  of  the  droplet 
was  uniform,  but  varying,  while  the  second  model  considered  tempera¬ 
ture  gradients  within  the  droplet.  Because  of  the  approximate 
nature  of  both  analyses,  and  their  rather  close  agreement  in 
predicting  heat-up  times,  the  bulk  of  the  comparisons  were  made 
with  the  infinite  conductivity  model. 


U8 

Regarding  the  heat-up  process,  it  was  noted  that,  at 
atmospheric  pressure,  ambient  temperature  exerts  the  primary  in¬ 
fluence  on  heat-up  rates  rather  than  preignition  chemical  kinetics 
as  is  the  case  in  bipropellant  ignition.  At  pressures  of  1  atm 
and  below,  the  data  seemed  to  substantiate  the  boiling  temperature 
ignition  assumption.  However,  above  1  atm  the  results  deviated 
significantly  from  the  theoretical  predictions. 

Further  study  will  be  required  to  determine  the  reasons 
for  the  failure  of  ti  .  heat-up  analysis  at  high  pressures  One 
possibility  is  that  the  final  wet  bulb  temperature  of  the  droplet 
may  not  be  well  represented  by  the  boiling  temperature  at  high 
pressures.  While  these  two  temperatures,  roughly  correspond 
at  low  pressures,  this  is  not  necessarily  the  case  at  higher 
pressures.  Also  it  is  possible  that,  at  the  increased  heat 
transfer  rates  at  higher  pressures,  the  temperature  profile  near 
the  droplet  surface  was  very  steep.  Consequently,  the  we;  bulb 
temperature  may  have  been  reached  at  the  surface,  while  the 
thermocouple,  located  between  the  center  and  the  surface,  indicated 
a  temperature  significantly  below  the  wet  bulb  temperature. 

In  the  combustion  study  it  was  found  that  the  experimental 

28 

burning  rates  did  not  correlate  well  with  Williams  theory. 

2  8 

The  major  assumptions  of  Williams'  theory  which  were  violated 
in  the  test  program  were: 


1.  The  activation  energy  was  not  sufficiently  large  to 
validate  the  thin  flame  assumption. 

2.  Natural  convection  may  have  been  important  in  the  case 
of  non-ad iaba tic  combustion. 

It  is  believed  that  natural  convection  caused  sfrioue  heat 
losses  from  the  flame  resulting  in  decreased  mass  transfer  rates 
from  the  droplet. 

Any  future  investigation  of  the  burning  rates  would  do  well  to 

28 

design  for  adiabatic  combustion.  Under  that  condition  Williams8 
analysis  could  be  checked  solely  with  regard  to  the  thin  flame 
assumption  for  a  range  of  activation  energies. 

Finally,  it  was  noted  that  the  ignition  delay  time  was 
significantly  longer  than  the  burning  lifetime.  In  particular,  at 
high  pressures  the  burning  lifetime  is  so  short  that,  after  ignition 
(presumably  at  its  wet  bulb  temperature),  the  droplet  is  consumed 
before  its  new  equilibrium  temperature  is  attained.  Consequently, 
the  entire  burning  process  was  transient  in  nature  and  the  mass 
transfer  rates  were  probably  lower  than  would  be  expected  in  quasi¬ 
steady  state  combustion. 


APPENDIX  A 


DERIVATION  OF  FINITE  CONDUCTIVITY  HEAT-UP  MODEL 


The  energy  equation  within  the  droplet  as  given  in 
Chapter  .1  is 
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where 


e  *  9  (sf  n)  . 


The  solution  of  the  ordinary  differential  Equation  A. 5  is 
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Differentiating  Equation  A.  7  with  respect  to  q 
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The  problem  now  becomes  one  of  obtaining  the  inverse  trans¬ 
form  of  Q. 
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:  This  his  a  simple  pole  at  S  ■  0 
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or 


Let  the  roots  of  this  equation  be 
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and  the  relevant  part  of  W*  is 
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APPENDIX  B 


B1  Fluid  Properties 


Property 


Density 


Specific  Kaar 


Latent  Heat  of 
Vaporization 


Critical  Pressure 


Critical  Tempera¬ 
ture 


Thermal  Conductiv¬ 
ity 


Heat  of  Combustion 


FLUID  AND  GAS  PROPERTIES 


Source  or  Method 

Value  of  Calculation 


86  LBM/FT3  at  20°C 

4.14  BTU/LBMR  at  20°C 

153  BTU/LBM 

35.8  ATM 

107 0°R 


FT  HR°R 

@  20°C 

-2000  BTU/LBM 


de  C.  Crater^, 
Beilstein^ 

Calc . -Johnson  and 
Huang^3  Method 

Calc.  Giacalone^* 
Method 

Calc.  Lydersen^* 
Method 

From  Vapor  Pressure 
Curve  at  Critical 
Pressure 

21 

Calc.  Weber 
Method 

From  Oxygen  Balance 

„  2 
Curve 
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B2  Gas  Properties 
Specific  Heat 

Heat-Up  Thermal 
Conductivity 

Burning  Thermal 
Conductivity 


~.56  BTU  at  2500°R 
LBM°R 


Cal.  Bennevitz  and 
21 

Rossner  Method 


Assume  k 


n„ 


air 


Keenan  and  Kay 


11 


.054  BTU  at  2500°R 
FT  HR°R 


Johnson- Huang 
Chart  using  Viscos¬ 
ity  Cal.  by  Bromley- 
21 

Wilke  Method 
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